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ABSTRACT 

The diabatic heating rate at 100mb is determined from the thermo- 
dynamic equation, utilizing previous calculations of the potential 
function xy and the "vertical velocity" w, for the period 1-5 April 1963. 
The four basic terms of the expanded thermodynamic equation; storage, 
stream-velocity, y-potential advection, and vertical heat transport are 
examined individually, and added to obtain Q, the diabatic heating rate 
which is compared to the results of previous investigations. The 
zonally-averaged heating rates are discussed in terms of stratosphere 


circulations and their associated physical mechanisms. 
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ive Introduction 

There has been recent interest in the diabatic heating rates and 
related energy changes in the upper troposphere and lower stratosphere, 
notably by Davis [1], Julian [2], Martin [4], K. Myakoda [3], and others. 
Diabatic heating is, of course, also at the heart of the problem of 
large-scale energetics. However, the usual approach through energetics 
yields average values for a hemisphere for the period of study, whereas 
the approach of this study takes the more direct route of obtaining 
12-hour average heating rates at each point of the usual octagonal 
(51 by 47) analysis-prediction grid system employed by the Joint Numer- 
ical Weather Prediction Unit, and other groups forecasting large-scale 
weather changes. 

Martin [4] demonstrated the feasibility of the grid-point approach 
by obtaining 12-hour mean values of both the potential function X of 
the wind and the "vertical velocity" G)= 2, Once these fields are 
known it is possible, utilizing the thermodynamic equation, to solve 
for on the diabatic heating rate. Martin's model is described in his 
paper [4] and only those aspects required for this paper will be re- 
viewed here. For more details regarding Martin's model, the reader is 
referred to his paper. 

The data period selected for the calculations was 1-5 April 1963. 
The data was provided by the Joint Numerical Weather Prediction Unit, 
Washington, D.C. This particular data period was selected in an effort 
to reflect the relatively large heating rates and the associated large 
meridional flows that normally occur during April, particularly in the 


polar stratosphere. 
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2% The diabatic equation 


If one expands the equation of thermodynamics, the result 


ae 5) Vv, xX | 

Sips oT + + Vy ‘VI+ Y, ‘V1 +138, = Q 
dt rae) ( (1) 

P P 
follows. Here Q is the heating rate per gram, while # and x give 
the nondivergent and irrotational parts of the wind in accordance with: 
VYekxvu+vx 
Here the fields of » and X were obtained preliminarily by Martin [4] 
by solution of the Balance Equation, and the xX-2 equation respectively 


(see section 3). When the equation of state and hydrostatic equation 


are employed in (1), we obtain our basic working equation: 


e(° se) + Wa REE) + Ve VEE) + ope Ye ‘o 


The fourth term of (2) may be transformed into the more convenient form: 


-TW=gUIdZ/|—-Ri 4s (22 
a se eee 


2P f (100 
= 4 sz 11240 + Lind (ZI) aN 


where 9 = ~1 39 but is also the negative of the coefficient of w 

in the final form of (3). “Here (2) and (52) are centered at /5mb 
APA AP 

and 125mb respectively, that is, at the upper-half and lower-half layer 

centers, respectively. The term -Q of (3) has a positive sign since 

both 0% /dp and the contents of the last bracket are negative. This 

agrees with the intuitive expectation that an updraft (w<o) is associated 


with vertical-motion cooling. 


The factor d2/3 mp is present in all terms on the left side of (2), 





and in order to determine its value at 100mb, the 50mb and 150mb data 
levels must be utilized. Unfortunately there was no vertically con- 
sistent 50mb data in punched-~card form immediately available, so it was 
necessary to derive the 50mb contour data completely from regression 
equations. In order to do this the United States Navy Weather Research 
Facility Regression equations of D. Lea [5] were used. These equations 


have the following form: 
= \7 
Bs = K + Ky (Z,,.) + a ue (4) 


where Ky, Kj, Ky represent the regression coefficients. It should be 
noted that the regression coefficients K,, Kj, Kg are listed by ten 
degree latitude bands from ten degrees north to the pole. It was 
necessary, therefore, to devise a computer program converting each 

Kn Cp) to a value of K, = Ky, (i,j) defined at each grid point within 
the 51 by 47 octagon. This conversion was based upon a subprogram 
giving sin in terms of the (i,j) grid. In addition K, CP) was made 
continuous with respect to latitude by averaging the values across the 
defining latitude boundaries for each ten degree band, while retaining 
the latitudeeband mean value. 

Since the 100mb contour analyses during the period of study were 
considered to be geographically consistent, whereas, those of tempera- 
ture were not, it was decided again to obtain TjqQ9 using T1900 = -£2e 
in (4). With this modification, the regression equation is mre P 


in the following form: 


Te / / / 
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\- Ky 





where 


/ 
Ky (Ky - 1) 


/ 
= Ke (gare 


(256 = 2100) = estimated thickness of the layer 100 to 50mb. 


and the term d2/3 np appearing in (2) has been expressed in the form: 


ae 
are ioe Deere) ar (Zioo7 Zico) 


Jinp Wms (6) 


in (5) and (6) (B50 _ Foo)” is estimated from the modified regression 
equation (5), with input data 2100 and (2100 - 2150) obtained from 
copies of data tapes provided by the Joint Numerical Weather Prediction 
Unit. The parameters #)99 and (#190 - #159) of (5) and (6) were used 
Wherever possible in order to utilize as much vertically consistent 
data as possible, thereby tending to minimize errors due to the use of 
the regression equations (4) or (5). All computations were performed 
using the CDC 1604 computer. 


Equation (3) can now be written in its final form: 
ae mere en <7 ae 
Cp 4] Bool + Lic ae C2 )T| Al O =e *Goa%e) 
(100 . 
Cp (2) wis Zeal asa} “Hie a He 6° be (7) 
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It should be noted that # = D + #, where saan is the standard atmosphere 


height and D = # - ae is the height anomaly relative to that of the 


standard atmosphere. In the first three terms of (7) 2, is independent 








of the variables x,y,t so that the time and space derivatives become 


alos Pash (Pax Ped]y Vf (Our Pee) * (Des Dis 


etc., and hence Zp can be omitted. However this cannot be done in the 
vertical motion term where , affects the vertical stability. 
Sic The irrotational wind and vertical motion structure. 

As pointed out in section 1, the development in this section follows 
Martin [4], who solves diagnostically for x for his  x-2 equation which 
is written: 
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oy (8) 
Th 


This equation is solved by relaxation using an iterative, scan-procedure, 


subject to the Bree oie. requirement: 


|x a < 4.0 «m 
Ky 


Then with (\ assumed to be known, one may compute 


A 
VV or ey x SP (9) 


Using w = w (p), as a cubic polynomial in p, with an upper boundary 





condition w = O at p = 0, Martin also derives an integrated form of (9) 


(100) 


giving w = w at each grid point. 


In terms of finite difference notation yw (100) has been shown [4] 


to be: (100) La 
Wore a Lb ah Vx i 
) =— J) (10) 


TNs 5 : : : : 
Here y is scaled in terms of centimeters, as this scaling was defined 





is (8). From (10), oo | takes on the finite difference form: 


(100) 


a 
S Sew 0% J 
Ud): : = ie V x (11) 


“) 


with Xi j assumed zero on the two outermost rows of the octagon. This 

was also the boundary condition used by Martin in solving (8). The 

47 by 51 octagonal grid mapping program used in the computations was furnished 
on tapes by the National Center for Atmospheric Research. Since it is 

not possible to compuceN Paton the "outer boundary" (Fig. 1), it was 


necessary to set w = o on the outer boundary throughout the computations. 


COMPUTATIONAL 
BounpARy 





Finite-differencing grid 


The normal five-point map grid 1-2-3-4-0 was used one row within the 
computational boundary, while a nine-point grid was used at all points 
ere inside. The nine-point grid simply inserts a three-point mean 
for a normal difference; for ae the aie i difference of a five- 


point grid is replaced with the weighted difference: 


(hare) +2 (es) + (A- He) /4 


in the nine-point grid. The actual five-point finite-differencing 





schemes used in the computations are as follows, relative to terms in ay /; 


(SE) = Sel k sty) - t- ; 
TU Sine)" aera (ts *)| By Sol 
(Y t, (5e)-(5 | 


VV Sin)” Te EY (Edel ms 
rete | (seae Ed, 


Furthermore the resultant Wij from (11) have been smoothed twice every- 
where interior to the two boundary rows just mentioned, giving ay the 
values used in this study. 

4. Results 

cai Be 


With the fields of Xij and known, as well as Vij and observed 
At “\ S 
from consecutive 12-hour map-fields, 0 has been determined from (7). 


Thestirst term of (7), callled the “storage term ae 


x || 
“eS care ( to0 a . (Zr. ~ To Sy (12) 


The second and third terms represents the 12-hour mean horizontal heat- 





transport at the gridpoint (i,j) and have the combined form: 





is TY [ested foo “TS +VAV aes He asa) (13) 


"py - transport" 'y - transport" 


The fourth term in (7) is the vertical heat transport (W)-transport) 


/ 028 2-25 ice aye TIZ46 + Lo al sea 


The four terms on the left side of (7), which have been listed in 


(14) 


the preceding paragraph, have also been computed individually as 
function$of (i,j) everywhere within the computational boundary. Com- 
bined values have been formed to give resultant octagonal fields of 
OG, )4 For a typical field of Q see Fig. 3. Zonal-averages of each 
of the five fields have been computed by programming an averaging 
routine for each five degree latitude band (over an area of 20-90 
degrees latitude). 

d, and the preceding four terms as well, have units of gpm which, 
for later comparison purposes must be changed to langleys. Q may be 
determined by a sum of terms of the form: 

Q= Mo. =o : fp 3 3d (9Az) _ M 3 3 3 42) | 
is 100 hes st (15) 
Here 0 is in ergs(gm)~! (sec)7! em? provided all terms on the right 


side of (15) are in c.g.s. units, Az=#~-Z in gp-cm and Mi 90 = mass of 








r] ° 


Fig. 3. Example of printout of Q values for (06 GMT, 1 April 1963) : 





the column 150mb to 50mb. 


: ee oe c 
= Me 5.419 X 10 ($e 0 (16) 
9PM ec 
The expression within the braces of (16) is now identical to (ST/st) 





in units of (°K) (sec). To get a corresponding value of 0 per gram 
per day, both sides of (16) are multiplied by 86,400 sec (day)! and 


divided by 102, the latter being the number of gms in 100mb: 


| 2 \ ad 
@(gen.) fe |s4i4 x16 X Sedco - 5PM, 


: 6st 3419X10 9,400 soo] A= 


0986 BPM coc 


Q = 6.298 x E (gpm) . eee) ly (gm)7! (day) + 


The units of 9 /3t [A#] are emphasized as gpm (sec)! because this is 
the form of the computer output units. Thus it follows from [3 (AB/e] 
=i 


that a geopotential change of 1 gpm (gm)~!sec is equivalent to a 


heating rate of: 
-| =] 
6.298 ly (gm) (day) 


The last statement is typical of each of the terms in (12), (13), (14). 
Zonal averages for each term and their zonal mean sum [0] of (7) have 

been determined and are listed by latitude belts in Tables 1 through 5. 
This is done also, for purposes of comparison, with computations of [0] 
made by Davis [1], who used a climatological or standardized atmosphere 


(Table 6). Note that the notation [ ] indicates a zonal average at a 


10 





specified latitude. 

In Table 6, the nine sets of results obtained here are averaged 
to give a mean daily average of the storage, the horizontal transport 
terms of both the stream and potential types, the vertical-heat 


1 


transport and the resultant Q, each in ly (gm) 71 day. Davis' results 


have been centered at 11/.5mb, but presumably this slight shift of 
levels is insignificant. In regard to the units of ly (mb) ~+ day! 
attributed to Davis, it should be recalled that 
= | cal 
|\O00 dynes 


whereas out unit of heat is: 


{ ly = | cal 
Se GBO dynes 


Thus the two heating rates, of this paper and of Davis, are very nearly 
compatible. The essential results are presented here as meridional 
cross-sections contained in Tables 1 through 6. 

In Table 6b, Davis’ computations for mean daily are based upon 
a climatological atmosphere at rest. It will be seen that Q should 
depend upon the nature of the meridional motion pattern, especially 
the field of vertical motion, in addition to any static distribution 
of radiating material. The term "Balance'’ is taken from Davis' paper 


and refers to balance requirements, on the basis that 
-(BR)= = Q = Co] Wy WY T+Vy: Yl-cu| 


with the "storage" very small (by several orders of magnitude) both in 


the present study and in Davis' 


Ee 





,-(w8) AT86Z2°9 = 1-235 wd3 [ :7-(Aep) ,-(w38) AT A0F AORVOeZ uoTsraauo0)] q-29s wd3 


[t-(Aep) 


OT J° 
sjtun ut *(¢€96T TtAdy G-1) autu ysnorzyi suo potiad AOZ qUOOT Ie sSejer aB3e104S paseroAe-A(Teu07 "IT oTqeL 


© iiinesel= yi lec Solel 6noeh nse°s 741 Vistar 
GSe°l- LOn°SI-  LOS*zi- zeora- 2e9°I ez*6l 9BB*LE | HLZ*LN | “Gdd=6 
| 
ee ee. NS... te SFT TEEERREET TET eR Perma re treeoees eee 
zn hEL*- 915° - Liesl 625 °E nooel | 
998°E- LlO°S= o6Brt-. Eales soars 60zteL SC eSBTLL oSETLL |b 
| 
a Pe ee e=tss, eae ee Be. 
eels SOL° ZE6°L bial zS0°l Lng°z | 
LiGe= ENS*E— Boz" HATZ LGLTS= —OzL ES” wots “atnt- fg 
ome em we oe en ee eo eS ee ——— oo ao oe oe oS 
S9S*l- 906°%-_—ELO°l-.—HOL S68°L 156° 1 d 
Bzhte-  ssse S96°l- d0c9- I162°- Seztol  9tztzzZkogtse | 2 
—_— = @@ as ae = oo => ——— = eee we we ee ee ee ee ee ww = oe i 9 8 oe we =e we ow we ow ee ee eww ee ew ee ee ew eee a ee ee ee oe oe oe oe ee eee 
99h" En GSL9° l= 6SL*L-  -Z80°l- —s_—«dBE°z- | 
nZZ°E OSE°E— LIG*A= = LZE*H- Ss EHOTOT= = PZSLE- SEOrLe— “ssetae-| 9 
We oe eei- ck 968° = teze 2 eee | 
ZnS*t-  «8h° Sind =L8°S 66S*tL = 6ONTLE—s 6 HBTON.—sAZES LN S 
ia Wi eeiiecso- | lliii0-z- ine*c- 902" =S=SssCGasetCtC(C(‘i‘C:ts 
gee 60E°9- sO ZOTLE- = BZL*BI- “SHS*E-  “E0L°9-  z2e0°0Z— 1EStzz-| 4 
— Dee tie? zcs- ss zce-t SS scoe> ~~ cuore. 
LzL°z EnZ°E g0L°- 3h0° 19n*6- . LOB'SI- 695°S olnrte | 
ac DOG deuen Tec, Se Sbe"c- | 75256-. 1 is 
n9L°h- =. 6S°9- OLN L genrE= nZara= tek tzt =" “L96*st= “ssere- | 2 
eS ee a oe ee oe a ee oe oe ee a ee ee ee ee eee ee oe. _..._.__-! _-__E! 
SSeEcL-_ sg LOL non nLLeZ £99°S 860°H | 
391°- SLS°2 yegrol —-usL70 = ASZ*E  “ozatle  “zssvol” clave | 
Hos-ss ss Gs-09—Ct—~«* 09-59  S9-02 Wisse SL-08 08-58. 8-06. 
ceed GZ-0€ O£=S¢ SE-Ih On-Sh 





12 





AT96Z°9 = {-225 wd3 JT :1030eF3 UOTSABZAUOD] 7_9eS wd3 ¢-OT JO S3yun ut *(€96LT Ttady ¢-T) 
auTU YySno1y} suo spotTaed AOJZ qMOOT Ie Sera Caiee Qeay AQTOOTeA-wWe9IZS pasexrsae-ATTeuOZ °Z7 2TqGeL 


—  gget-_ — Z08 THE —sCOOAHTZE SEH OZ ~—szaL°O9-.s-ghl°So.  SqVLOL 
Lnn°Bz =«-PES*BLE «—SeETIEH- Loz*ON- “G9Z*LIC- Searshi- legen? “OSo"LZE | “dud-6 
oe ee ew oe =_ ke Re a a oe ee ee a a 

699°9 Z51°9 h6H*OL-  _zSO°E SOL*L Lens [ 
nSS?l ZSS°S9° | LOB*HDL- “SEOTLO SHEcZTH- OSL*EE- Seervl2 DIETLS | 6 
nip ee ae ea ie RPE ee a 2 
Hh eh 132° 968°S- S92°L 969°€- _—hbO°ES r 
ELL°LL «© ZELTLL = BS9TOLI- B95°ES —LOD"SI- ‘sansi- O99°9Z BnzrsL | 8 
| 
a ee a a eee ne: 
pecs 99E°L \eb2e— ie ice 959°z-  ___300°8- | 
pateeecOG  neze6=  “Gogest Ji6°oe- 966°- ezstl=zngrS S| 
— =e =e om oe cee coe oe oe ee we wee ewe ee ee ee es es es es ee es oo ee ee ee es ws ws ws ws es wes we es es nes es ss we es es ee ee ne es es i es es ee es is es se es is we l — ome oe oe oe oe o> 
OnE °— 91E° Lon? t BSE°S 938°n- BHO i 
Sen*LL = SLL*BL = PON*LE=-LSOT9L-  EBOTLE- EETEL- OSz*L- ‘LnBSE | 9 
| 
eee cs el 
HZ6°h- _ZBLER 189°SL RHEE Elnes- ZL LL | 
961°9- 968°9G 9LOTSE- L90°HH- “B9D°ZL-. SBS*LZ—- “O9L°H Ise"so |G 
I 
eee ape ee ee. Mk dl 
982 °2 Lng° FROTLL ENT 966°2L- 32S°ZI ! 
HEO°SL- 6£9°2 acces 619°9 Gen°9E- LO0°9@- CZh*he  -AAzZ*BS | 
—_— = = == == ee ee ee = GE GR CEP cm oe ce ee ee Oe ee ee Se ee ee we ee i a oe Ee eee = eee ee, me ee se ee oe oe 
ZS9°L 906°2 810°S n90°E Z6O°LL- HHS 7 
hO6"9- LOOTEL  SEneHL- OhHeEL- “HSI*SL- ‘BSSe*e- LLS*6 mine9z | & 
— ce ce es a es a ee es es es es ne es me GR ce ww ws es es es es es os es i es i es ws ws ws es ses ss es ee se i ‘ ee ee es 
sep ole |e eSh BEL °h B90°E Lig°2- hoe DL . 
E60°l- « EBO*E- —HOVTOE «= -S9TLA~ ZEBTOL «= ELETOL= L60%L- — LontHL |Z 
Rb ec pool SE <b Dn gr 0 ge eee 
eezae mennolc 1 n68°S aL) i ee ee [ 
PSO ISh?Ofs | GGs=e-9  PessOh- OF6°9S- L¢s°6- 2¢s°z- 
bees Seat) NLS) SS) a 7 ee Siu Me Gace | EC Ean © 


OC=S2 Gc-0f O€—-SE S&-On On-Gh 





13 








[y-Aep 7-w3 41g6z°9 = 7_9es wd3 | 


:10}208J UOTSIBAUOD] 7_o9esS wds c-O1 JO s3qtun ut 


‘(€961 Ttady ¢-1) ouTU Yy8noiyA sUO spoTied AOF quQOT 32e UOTROeApe [eTIUs}Od-X paBerzsae-ATTeUOZ “°F JTqGeI 
—— eit, 898 HL=, gaz°S|> | [p6°e=,_, 2gl°9e, | azl°09=_, | S¥wiOL 
OO6°LB- RSL°OH- LOZ°SB- SEO°hhtL SLH°982 E2S°6C2 Ben°sdt z29s°es | Aud-6 

1 
ee a ee ee ee ee ee ee ee SS ee ee eS eee ee ee -|-------- 
Gic° Z60°9L- OL °9L-. dL 2zP°n- SSE °L- Zl0°Sl- | 
G2i°Zl OFL°Sl SL2°S9 3° oH BEL ac9°cl nc6°Sl- LEscl- 96L°c- | 6 
shoes. ------____---- : ee 
E91 °h- OZ6°Eh- SiB°9 ELge- ne 388° | 
O9L°LH- fl49°OL- O16 %HH ZO9°E Fa0°l ZHO°EL OL0°21 Lane | 8 
ee Sa Se ee a ee ee eee ee Se a ee a ee ee Oe es 
gSl°l OFL°L- _ SOM*LZ- ZLLeZ-— OE Faust | 
G6S°h- BOE*ZLI- OL *LLI- SLz°9Z2 661°69 nBL°9Z SEG BORE y] 
SS SS SS ES SS SS eS a ae ee ee a ae a es a ee 
PEO 5$90°86 866°LS Seco OLL°- 5LO°EL- ! 
6E2°2 £96°06 OSL2°L2- = LBA L- 5382°89 Er? OL hEO°LZ Sane Ol | 9 
a es aoe ee ee ee ee ae ee eee ee eee CC 
162°S LSh°0E BRL *- heneL- OL ° 782° nL - fi 
68L°SC- BOL *HH- ZEE*SH- FOR LI L9h°cl LS9°¢s heO°sz BESS St G 
_— oe oe ew oe ao oe Co ee eo oS So ee a eee aaa ee a i ee os 
z50°2 So9°LL = SHL*L Lint heh°S-  F66°HL- , 
ZLB ° He S$82°1L6 L6B°09 998°92 L9E°6S FLL °69 B8L°9S gene t 

— 77 """Nenege. o8e°Zhe. Len*LZ=.269°=. 082" ne e2e 
CBT°Se- = OK OD n9O%HL Z9S°6¢ HEE eS Z60°SS- 992°8Z- 6LS*2L- | € 
5 a ei Sa iy a ee ee ee ee a a a ee ee a ee ee ee ee so Se = 

Seerl- LEZ S60°lE —SBL*h asi°l 2n9°9 , 

BLE°ZL- BSL°08- EOR *HH- LZETOH O1e°SS OLn°Eh FZE°9 | 

ni ee ee ee eee ew we ws ee 
n9L°- 615°6S-  €02°LH-  _30L°- 6LE°I : . 

ehh? L- J62°ESG- ZBlL°bS- Llarzl-  nSGl°z 9£9°02 00S°22 
! 

‘0S-sS SS-09 09-59 S9-o2 Diss 5S 51508 PSEA  SAaaa 

02-SZ GZ-O£ O£-SE SE-Of On-Sh 





14 





[,-AeP tous A1862°9 = {-225 wd3 [ :1039eF UOTSIZAUOD] {-225 ud3 t-OT JO satun ut 
“(€961 Ttady ¢-1) eutu y8noiy} suo spotisd 20F quooTt ye ArAodsuerz, Jeoy [PITIAVA posersae-A[[eUuoZ ‘*» 2TqP] 


~ [ee"e_S8176- | g9°E1 =, 6tS*1=-, 18L79, | -225°- rstvio1 
601° ZS8°- €09°Z- 88h hlL- k69°E- 798°SBl 9L6°9r 820°2S | Gud-6 
| 
mee eee ee — ene ee ee eee ee eee em] ww mw ew eo ew ew ew ew ew ew www @ @ ww ww ow = om ow = we = oe we ww = ow www oo ww o> oo oe = ow oe 
GLL°ZL- = LOL RR LL9°- 516°- OfF0°Z ahi? t 
USA S6L°- 9EL°- enSG° eyLss Sol} ZLEcL- SEG oa= St 6 
| 
a SS Sr SS SS iS SS SS SS SSS SS ESS SS 
8SL9°ZL- —s--O 9 L O° hh 88z°9 OZ0°z- neo a Bae 
I9L°- LiS°- GLL°- O6t° | HLse- BOL°I- OS6°I L9z°9 | g 
| 
cee ae ase Se ee IS SS eS SS SSS SS SS = SS SE Se 
716° LeusGe EZ0°L- GGG Psc= [on°- 
69E°%=- 8LO° SEL°C—- 62E°9- ciz° Li8°2 6Sh°Z L1¢8°S- | 1 
| 
ye eo Se Se SS Sa SS SS Se SS SS eS SS EE Oe Se ee ee ee eee 
190°S2 982° 60h°ZL- OS6°L- Cie Lin 
56E° LOn° 698° 819°Z- 696°C- 98h t 66hn° Ol 6ZS° ll | 9 
een ee = ee eee SS 
LOZ°E hth? Z- Z6E°L- 682° 9n0°Z 980°L- 
SL0° L9Z°- ZLB°- O16°L- hOZ*L—- ZL2°- hl 9°S 098°6 | G 
| 
SS SR OS OS CO SSS SS OO OO Oo SS OS SoS Ss SB Oa SS eS BS Se ee a ee a a ee ies SOO ao 
. EHO°OL 99n° hEQ°h-— L9O0°E- 162° 929° L 
08s? LSt° 9GzZ°l LELoL- G90°C- 56h el Lhhect GOS*9Z | t 
| 
———— ee een 
LLS*tl- LZl°Ee- esorl hoses 91L0°- JLO°L- 
L9¢° 6Sh°- 062° BOL°L hnicZ°e O6L°Z 6Sn°S- B80E°Zl- | ¢ 
| 
Se SO oS SS ee we we we ww ww we ww ww es we we we we we we we wm ae a ae i a ew we ee we ee wm we ee me es i ee we ee mm we wm of ee ee ee ee ee we 
NOLS 6L1°9 LOO°L 97° h- 160°L- 960°- 
ZBL°- 189°- 8Z0°?- hl Z°e- FLZ°- OL1S°h GLZ°9 869° Z 
Bs oe eee eS fs ere 
LZL°Bl- =: 9:09 °6- 62S°S O9E°S ihZ°Z LEz° 
Lhnz°- LEL° C8C°- 806°C- LZ3°Ee- ZLO°C- B80S°h GAOS HL | l 
| 
0S-SS GS-09 09-S9 S9-O0L OL-S2 SL1-0 98-S3 G8-06 ud 


Oc-S¢e 3 =e O£-S& S¢-Oh On-Sh St-OS LVI 


15 





[,-u3 t-AeP AT86Z°9 = [~99S wd3 [T +:10}0eF UoTSIPAUOD] T[_99S wd3 701 JO sjtun ut 


“(€961T Ttady S-1T) seuTU Ysnory4) 9UO potied 10J qUOOT 3e SsazeA BuTJeDY DTReQeTpP pos8er9aae-ATTeu0Z °C eTqeL 
ee ee ee oe a ae eo ee le a ee a ee a a a oe a ee ee eee ewe ee eee 
9S¢°E- 9nh?S— ESE°ElL-  SLE*L- 6ZL°9 180° 1 siviol 
566° eE5°C FLil°9- iS6°rti- NS ee rOn° lt COL le B8C°SS | OYd-6 
ae pe ee eee ee se -__ a ee 
BHO°?L- r69°h- FEL ch6°- Lho°? rSt°l 
9% ° gon? 98L*L- F9OLPL OFS °¢ FLe°l CEO. 1 — Coy. i= | 6 
= yeh ee ee a ae 
6c9°Cl- RSO° fh OSC°9 200 7c £00°= 7 a 
S00Q° 90c° Eco l= O£0°C ACG | oes 9G. la OC"? Fines | 8 
ae ee en, a a ee ee ois bo ek 
016°8 h9O9°S- FSO°L- RLSes ck ° tS? = 
6¢0°- BcS° cS$0°*e- iSL°9- ESO; — 298°2 R8zx °C 008°S- yh 
i 
_ SORE a ans xO ee Ee Se OS SESE Se ae eee ae a een 
h90°SC S8c°hR S9E°cl- 968°l- £6c° 5ith° 
hls ° Bell LSo°- GShB° Cc Ghte °S— rSe°rR 98n°Ol LB38°ILt ! 9 
i I 
ee ee ee Se a ee ee eee ae ee oe oe oe ee a ae eee a 
CUS Ese C— Gec°l- gece ° S55" 1 §90°1L- 
100° 380:° Coit" t= fins c hen°?t—- B88h°- 799°S O0cS°Ol ! S 
| 
= ee aw © aw oP oe ow om aw ow oe ow ow ow 6 Ow ow ow oe ow ow ow oO ow © ow © Ow ww © ww ww ee ww = ow © ow ow ce ce cm co ce ce oe Om oe ewe wer oe a @® ow oo of @® ao ow oo @ ow © @ @ ww ww & ow fr wrwrwrwroen-— 
999 °91 son? Ll hlo*h- GcO°s- Zi9° rai cas | 
Ogh® - WRBE® 6S0¢°t n9O°L- Veh’ c- rzxco°l 98l°cl 280°12¢ | th 
| 
— = aw ow ao ow = a _ om cm ce cm ee oe ew Go aw ow aw om Oe oo ow oO ew cep 6 oe om cm cm oe SO ow ow ow ow ow ep om om ow om ow om ow ow ow ow ow aw oe oe —am ow em om ae we ow ww a ww i 8 2 2 2 Pe = 
09S °hLl- Cowes — he lel GSS°S JO, gee Le Sicl= 
86c° bce °— 99° R1S°L Z0L°¢S h60°Z e9E°S- hhrQO°cl- ¢ 
' 
= = aw ow aw om eo of oe om ow cw ow Ow ow ow ow ow ae ee See ew 0S SS S22 222 P22 PP 29 OF P22 2 O22 2 2 Oe 2 TFT 222 22 ee 2 ee Teo 
O6L° I51L°9 _B8hOL EEceh— L2L°l- ILO°- ! 
F6t°- 3LL°- LCL A= SSC oe hOl°*- Lon? t hhl°9 CRB°s | c 
Se i St SS SS = SS SS a a a ae a aa = a a 
S61C- Gus E£eS°S- BL9O°ES che°’s G9L°C 5SS° 
9thke- cOS° 9¢2°- 9LE°S- OTA a= 891° c- O8h? t n9O°St ! L 
| 
OS-SS SS-09 09-S9 G9-02 Ol-S2l ¢1-08 08-S3 


; 02-SZ GZ-OF O€-GE S€-Oh On-St 





1G 





45,|45-65 [50-56 [55-60 | Go-BI6S-10 | 70-75 | 75-90 | 80-85 |85-Fo 


| BALANCE. RQ] SD 


Q 


a 
do 


4 


17 


THERA LED 









CooL ute 
BALANCE REQ], OG | .08 





13 


Uf 
“~1) 


OF 
4 





See 


ant 


4 bs [e16 Is17 Fa 14 


(b) 


24-hour averages of storage, 


7 
& 


=n 


Q Bow amor 


advection, xy-advection, vertical heat transport 


p- 


(a) 


Tadeo 


i averaged meridionally 
Q of this study with those of Davis (ele 


ions to the daily .diabatic heating, 


and their contribut 
Comparison of 


(b) 





aye Conclusions: stratospheric circulations and associated physical 
mechanisms 

As already noted, the zonally-averaged storage term is very small. 
In addition the stream- and potential-advections, averaged over five- 
degree latitude bands are smaller by approximately two orders of magni- 


tude than [Q]. Thus to a good approximation, valid everywhere at 


least as to sign, Qle —¢,(7 U9) 


and is of the proper order of magnitude as well. 


By the equation: 


loo EEL I=-815V wp 


it follows that [w 99] >O occurs when [W.\y] represents convergence. 
Now convergence occurs near the pole, since the expresston forstnc 


zonally-averaged divergence 


[Vo] =[S4 — EF [Se] tang hn 
M44 


(18) 


with the second term (which represents the convergence of meridians) 

becoming dominant [Martin, 4] at latitudes > 75N. At lower latitudes 
[VV] = Xieo| = Sol /33 

is a valid approximation; that is, the spherical divergence is only 

slightly affected by convergence of meridians. Hence, where F100) is 

negative, there is divergence, if boy oof is positive, there is convergence. 


Consequently at 100mb, we have centers of convergence at the pole and 
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latitude 35, and of divergence at latitude 65. 

In inferring a stratospheric circulation, one may start with an 
interhemispheric meridional circulation proposed by Kellogg [6], 
especially applicable, in some of its aspects, to the winter-summer 
reversal in the stratosphere. The octagonal grid, together with the 
associated boundary conditions, is not able to resolve the effect of 
the cross-equatorial flow to at least latitude 20N latitude. However, 
because of the decrease of mass in the winter stratosphere, cross- 
equatorial flow must occur near the time of the vernal equinox. This 
evidently is responsible for the convergence center near latitude 35N. 

Likewise with the return of the sun to the Northern Hemisphere, 
ozone regeneration goes on rapidly, followed by ozone-heating, 
especially at 45 km. According to Pressman [7]. an ozone radiative 
heating maximum proceeds progressively from 30N latitude to the pole 
between the dates 15 March and 15 May. For time of this study, 1-5 
April, 1963, one might expect maximum ozone heating at 45 km at or near 
60N latitude. This effect is so marked as to cause a strong forced 
outdraft (Cy YO) at 60N latitude, and at an elevation of 45 km. In 
order that three-dimensional mass divergence be within reasonable 
bounds, there must be a sympathetic updraft from below. This is the [w] 
measured at the l100mb level and presumably exists at levels between 
100mb and 45km. Also according to Martin's model, sizeable w-values of 
the same sign still persist at the 200mb level, but are smaller than 
W199 in units of cm sec~!. The level of "return flow" in these cells 


cannot be designated from presently known studies. 


ry 








Fig. 4. Proposed meridional flow patterns 
early in the spring~-reversal period 


The "maximum" convergence in the vertical and associated downdrafts at 
the pole at this particular time of ozone activity should then be 
viewed as a "spill-over" effect from 65N latitude. The maximum con- 
vergence at or near latitude 35 North, which mirrors a zonally-averaged 
100-mb maximum values of &)at this same latitude is postulated as the 
primary result of cross-equatorial flow. Considering the level at which 
the greatest polar-vortex intensification in the Southern Hemisphere 
occurs, the level of the 35N latitude maximum of convergence has also 
been placed at 45km. 

Fig. 4 together with the above discussion "explains" the mean 
meridional flow pattern in the stratosphere during the period of this 


study. One may offer as conjecture that as the stratospheric summer 
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advances (say by June 15), the maximum convergence will have advanced 
to 45N latitude while the zone of maximum divergence associated with 
most effective ozone-heating near 45km advances to the pole, there to 
become associated with divergence at higher elevations. 

As contrasted with the 0 of this study, Davis’ results show only 
moderate cooling at every latitude due to predominance of long-wave 
radiation over heat sources. Two features by which April, 1963 
appears to differ from climatology are shown in Figs. 5 and 6, both 
taken from analyses of the Free University of Berlin [8]. The latter 
figure shows a pre-existing zone of vortices at latitude 60 North, 
as of 1 April, 1963, 1200 GMT. The other figure mentioned above shows 
that, as compared with a completely different April (1962), pressure 
was higher over the pole but lower than in 1962 in a vortex ring at 
60N latitude. Fig. 6 shows a possible reason for this: strong baro- 


clinicity on the equatorward side of the vortices, with ‘zonal 





————_—— TT 


+ (6%) [al 


computed over 45-70 degrees North latitude favoring the continuation of 


generation’ 


such a baroclinic zone and westerly winds. In fact, according to the 
analyses of the Free University of Berlin, an easterly regime with an 
associated polar anticyclone did not occur until May 15, 1963. 


Finally, how may one reconcile the zonally-averaged equation: 
Q-(- edu) arn ® 


so that both horizontal flows and storage are relatively unimportant 


C19) 


ZA 





at Samm. 
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Fig. 5. 100mb height changes, April 1963 - April 1962 from analyses of the 
Free University of Berlin ae 
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as determined here. Equation (19) implies that with —6{Td}>0 (down- 
draft motions) , |Q]<o is the result. Such a result is reasonable on 
physical grounds by noting that descending air is warmed dry adiabati- 
cally, and in addition one may expect some increased diabatic solar 
heating acting upon an ozone increase in lower stratospheric elevations. 
However, one immediate opposing result is that the C05 in the warmed 

air will emit at a greater rate, tending to re-establish the heat 


balance. The same argument applies in reverse for updraft motions. 
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